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Abstract 

A  model-based  LQR  method  for  controlling  vibrations  in  cylindrical  shells  is  presented. 
Surface-mounted  piezoceramic  patches  are  employed  as  actuators  which  leads  to  unbounded 
control  input  operators.  Modified  Donnell-Mushtari  shell  equations  incorporating  strong  or 
Kelvin- Voigt  damping  are  used  to  model  the  system.  The  model  is  then  abstractly  formulated 
in  terms  of  sesquilinear  forms.  This  provides  a  framework  amenable  for  proving  model  well- 
posedness  and  convergence  of  LQR  gains  using  analytic  semigroup  results  combined  with 
LQR  theory  for  unbounded  input  operators.  Finally,  numerical  examples  demonstrating  the 
effectiveness  of  the  method  are  presented. 
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1  Introduction 


The  use  of  shell  models  to  describe  structural  dynamics  is  pervasive  in  applications  ranging 
from  noise  reduction  in  aircraft  to  flow  control  in  flexible  pipes.  While  general  shell  equations 
can  be  used  in  a  variety  of  geometries,  they  all  share  the  property  that  component  displace¬ 
ments  are  coupled  due  to  the  geometry.  This  leads  to  significant  challenges  when  developing 
appropriate  models  and  approximation  techniques,  and  constructing  effective  controllers. 

In  this  paper,  we  consider  cylindrical  shells  due  to  their  prevalence  in  applications.  Con¬ 
trol  is  provided  by  piezoceramic  patches  bonded  in  pairs  to  the  surface  of  the  shell.  These 
transducers  provide  significant  actuating  capabilities  due  to  the  piezoelectric  effect  in  which 
input  voltages  generate  strains  in  the  patches.  Utilization  of  the  converse  piezoelectric  effect 
(strains  produce  voltages)  also  allows  the  patches  to  be  employed  as  sensors.  When  combined 
with  their  light  weight,  space  efficiency  and  reasonable  cost,  these  properties  make  the  patches 
highly  effective  control  elements  in  a  variety  of  applications.  From  a  mathematical  perspective, 
the  use  of  surface-mounted  piezoceramic  patches  leads  to  unbounded  control  input  operators. 

Experimental  work  has  already  demonstrated  the  potential  for  success  when  employing  the 
patches  as  actuators  in  applications  involving  cylindrical  shells  [8,  13].  However,  these  initial 
investigations  have  not,  in  general,  utilized  the  full  potential  of  the  patches  due  to  limitations 
in  hardware,  models,  approximation  techniques  and  control  laws.  For  example,  a  common 
means  of  calculating  control  gains  is  through  the  use  of  modal  expansions  [10].  However, 
closed  form  expressions  for  the  modes  can  be  determined  only  for  a  limited  set  of  models  with 
severely  restrictive  boundary  conditions.  The  use  of  incorrect  modes  when  calculating  control 
gains  can  lead  to  loss  of  control  authority  and  possible  controller  instabilities.  Hence  for  most 
physical  shell  models,  modes  must  first  be  accurately  approximated  if  modal  methods  are 
employed  for  control  design. 

In  this  paper,  we  present  a  model-based  method  for  controlling  shell  vibrations.  For  sim¬ 
plicity,  the  Donnell-Mushtari  shells  equations  with  Kelvin- Voigt  damping  are  used  as  a  model 
(the  assumption  of  strong  or  Kelvin- Voigt  damping  is  reasonable  and  typical  for  many  shell 
materials  such  as  aluminum).  The  methods  are  general,  however,  and  can  be  applied  to 
higher-order  models  (e.g.,  Byrne-Fliigge-Lur’ye  model)  if  the  application  warrants.  A  general 
Galerkin  method  based  on  splines  is  then  used  to  discretize  the  infinite  dimensional  system  (see 
[7]  for  details  regarding  the  numerical  method  and  a  comparison  with  finite  element  methods 
for  shells).  Through  the  choice  of  basis,  the  method  is  constructed  to  be  flexible  with  regard  to 
the  boundary  conditions  and  material  nonhomogeneities  which  arise  in  typical  applications. 
Furthermore,  development  of  the  model  and  approximation  method  in  terms  of  a  weak  or 
energy  formulation  facilitates  consideration  of  the  distributional  derivatives  which  arise  when 
including  patch  contributions  in  the  model.  This  provides  a  setting  suitable  for  direct  simu¬ 
lations  and  control  design  as  well  as  computation  of  frequencies  and  modes  for  the  shell. 

The  model  and  approximate  system  are  then  employed  in  an  LQR  full  state  feedback 
theory  to  obtain  feedback  gains  and,  ultimately,  controlling  voltages  to  the  patches.  While  full 
state  measurements  are  not  available  using  current  instrumentation,  and  hence  the  techniques 
cannot  directly  be  implemented  in  experiments,  they  provide  an  important  first  step  in  the 
design  of  effective  compensators  based  on  state  estimates  calculated  using  a  limited  number  of 
observations  (see  [5]).  The  consideration  of  the  LQR  performance  also  illustrates  properties  of 
the  system  and  model-based  control  techniques  and  facilitates  investigations  regarding  issues 
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such  as  patch  number  and  configuration.  Finally,  the  consideration  of  the  problem  provides  a 
step  toward  the  development  of  model-based  controllers  for  fully  coupled  structural  acoustic 
and  fluid/structure  systems  involving  cylindrical  shells. 

The  strong  and  weak  forms  of  the  Donnell-Mushtari  shells  equations  are  outlined  in  Sec¬ 
tion  2.  In  presenting  this  model,  care  is  taken  to  include  both  passive  (material)  and  active 
(actuator)  contributions  due  to  the  patches.  An  abstract  form  of  the  model,  based  on  sesquilin- 
ear  forms,  is  also  presented.  This  provides  a  natural  setting  to  prove  model  well-posedness  and 
convergence  properties  of  the  LQR  control  law.  LQR  full  state  feedback  laws  for  systems  with 
no  exogenous  force  or  forces  which  are  periodic  in  time  are  presented  in  Section  3.  In  the  former 
case,  convergence  of  the  approximate  suboptimal  gains  to  the  optimal  gains  for  the  infinite  di¬ 
mensional  system  is  proven  using  analytic  semigroup  theory  in  combination  with  LQR  results 
for  unbounded  control  input  operators.  A  Fourier-Galerkin  method  for  approximating  the  sys¬ 
tem  dynamics  is  outlined  in  Section  4,  and  the  effectiveness  of  the  LQR  method  for  periodic 
forces  is  demonstrated  through  a  numerical  example  in  Section  5.  This  example  demonstrates 
that  through  the  use  of  the  model-based  methodology  with  general  Galerkin  approximations, 
significant  attenuation  in  shell  vibrations  can  be  obtained  using  piezoceramic  patches. 


2  PDE  Model 

The  system  under  consideration  consists  of  a  thin  cylindrical  shell  with  surface-mounted  piezo¬ 
ceramic  patches.  It  is  assumed  that  the  patches  are  mounted  in  pairs  with  edges  aligned  with 
the  circumferential  and  longitudinal  axes  of  the  shell.  The  edges  of  the  shell  are  taken  to  be 
fixed  in  accordance  with  common  experimental  clamping  techniques. 

To  specify  the  geometry  for  the  corresponding  model,  we  consider  the  longitudinal  direction 
to  be  aligned  along  the  x-axis  as  depicted  in  Figure  1.  The  displacements  of  the  middle 
surface  in  the  longitudinal,  circumferential  and  transverse  directions  are  denoted  by  u,  v  and 
w,  respectively  while  the  length,  thickness  and  radius  of  the  shell  are  denoted  by  £,  h,  R.  The 
region  occupied  by  the  middle  surface  is  denoted  by  Fq.  Finally,  the  shell  is  assumed  to  have 
mass  density  p,  Young’s  modulus  E,  Poisson  ratio  u,  Kelvin- Voigt  damping  coefficient  cc  and 
air  damping  coefficient  p. 


x 

Figure  1.  Thin  cylindrical  shell  with  surface  mounted  piezoceramic  patches. 
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Actuator  and/or  sensor  capabilities  are  provided  by  s  pairs  of  surface-mounted  piezoce¬ 
ramic  patches.  It  is  assumed  that  all  the  patches  have  thickness  hpe.  Young’s  modulus  Ep^, 
Poisson  ratio  i^pe  and  Kelvin- Voigt  damping  coefficient  Furthermore,  it  is  assumed  that 
the  glue  bonding  layer  provides  negligible  contribution  to  the  structural  dynamics  (the  reader 
is  referred  to  [3,  5]  for  details  concerning  the  incorporation  of  differing  patch  characteristics 
and  bonding  layers  in  the  ensuing  models).  The  region  covered  by  the  patch  pair,  with 
edges  at  Xii,X2i,  is  delineated  by  the  characteristic  function 


Xpei{x,0) 


1  ,  Xu  <X  <  X2i  ,  9u  <  9  <  92i 
0  ,  otherwise  . 


The  indicator  function  Spei{x,9)  =  Si^2{x)Si^2{9)  ,  where 


1  ,  iC  < 

0  ,  X  = 

—l,x> 


{xu  -f-  X2i)l2 
(ajii  -F  X2i)l2 

i^xu  "h  ®2«)/2 


1  1  9  <  {9u  +  92i)l2 

0  ,  9  =  {9ii  92i)f2  , 

— 1  ,  9  >  {9ii  +  92i)f2 


delineates  the  sense  of  the  forces  generated  by  the  pair.  The  symmetry  of  the  function 
arises  from  the  property  that  for  homogeneous  patches  having  uniform  thickness,  equal  but 
opposite  strains  are  generated  about  the  point  (xi,  9i^  =  -f  X2i)(2,  {9ii  +  ^2i)/2). 


2.1  Strong  Form  of  the  Modeling  Equation 

We  consider  here  the  modified  Donnell-Mushtari  equations 


dt^ 


dx 


89 


2  =  1 


P  ,d^v  dNe 


dt^  89 


8x 


Rn  ^{^e)pei  c  (  a\ 

RQO  q0  Rpei{Xj9) 


i=l 


8‘^w  8w 

Rph-g^  +  Rp-^ 


R 


8^M^  1  8‘^Me 


8x^  R  89^ 


-2^  +  Y. 
8x89  ^ 


Rqn-Y^ 


i=l  L 


R 


8^{M,)pe,  ,  l8\Me) 


8x'^ 


+ 


pei 


R  89^ 


(2.1) 


as  a  model  for  the  thin  shell  dynamics.  As  detailed  in  [3,  5,  9],  these  equations  are  obtained 
through  force  and  moment  balancing  with  only  low  order  terms  retained.  Here  Mx,Me,Mg^ 
and  MxB  are  internal  moments  while  Nx.,Ng,  Ng^  and  N^g  denote  internal  force  resultants. 
External  surface  forces  are  denoted  by  qg,  qn  whereas  the  external  resultants  (line  moments 
and  forces)  generated  by  the  patch  pair  are  designated  by  {Mx)pe,,{Mg)pe^,{Nx)pei,iNg)pe^. 

Expressions  for  the  internal  force  and  moment  resultants  are  derived  under  the  assumption 
that  stress  is  proportional  to  a  linear  combination  of  strain  and  strain  rate.  This  yields  a 
model  which  incorporates  Kelvin- Voigt  or  strong  internal  damping.  As  detailed  in  [3,  5],  the 
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resultants  N^;,  N^e-,  Ngx,  M^e:  Mqx  derived  under  this  assumption  are 
Nx  =  — 2  X]  TTZ — 2~  Xpeii^t 

^  ~  ^  i=l  ^pe 


coh 


pe 

2cx)  hr 


_  ^  (4  +  ^eo)  +  Y,  -.  Y  f 

1  2=1  ^pe 


N,b  = 


AT-  I  Ep^hpe  . 

Nqx  —  TTTi  :  \^x6  +  2^  Y^  ^^xeXpeiy^-i^j 

i-\  v-^ 


2(1  +  u) 


+  ^pe) 


Cvh  ^  I  ^Dpe^pe  ^ 


-£x9  +  X! 


^SxOXpSii^'! 


Mx  = 


2il+u)  ^(1+^/pe) 

Eh^  ,  ,  ^  2Epea3  ,  X  /  /IX 

(Kx  +  V'Kff)  +  ^  (Kx  +  ^pel^e)  Xpei  {^i  ^) 


12(1  -  Z.2)  ^  ^  ^  3(1  -  ulJ 

+  12(T'- y^)  '*'  '*'  ^  ^  Xp.,  (i, «) 


(2.2) 


7!  /T  7i  /T  Eh  Ep^Cl^  /  «x- 


24(1  +z/)  ^^3(l+7/pe) 

where  the  constant  a^  =  (/i/2  +  /ipe)^  —  results  from  integration  through  the  thickness  of 
the  patch.  Expressions  for  the  resultants  Nq  and  Me  can  be  obtained  by  replacing  Sx,  £$■,  i^xi  >^e 
in  the  expressions  for  Nx  and  Mx  by  ee-,€xiKB,Kx,  respectively.  The  midsurface  strains  and 
changes  in  curvature  for  the  Donnell-Mushtari  model  are 

du 

F)t.  n  (in  n  (i:i:  n.  (in 

(2.3) 

d^w  1  d'^w  2  d'^w 

Ke 


\  dv  w 


dv  1  du 

^  ^  'Rdd 


Kx  — 


dx^  ’  ■■'''  dO^  ’  Rdxde 

Note  that  for  the  undamped  shell  which  is  devoid  of  patches,  the  resultant  equations  (2.2) 
reduce  to  the  classical  Donnell-Mushtari  expressions 


Nx 


Ne  = 


Eh 


(l-i.2) 

Eh 

(1  -  1/2) 


Nxe  =  Nox  = 


du  u  ( dv  '' 
\  dv  w  du 


Eh 


Mx  = 


Me  = 


d^w  1/  d^w 

+ 


dv  1  du 

2{l+u)  [^^R'^\ 


,  Mx0  =  Mbx  = 


-Eh^ 

12(1-1/2) 

-Eh^ 

12(1  -  i/2)  [i?2  06^  '  "  dx^\ 
—Eh^  d'^w 


dx^  R?  dd^ 

1  d'^w  d'^w 
^v- 


12.R(1  -b  i/)  dxdB 


(e.g.,  see  [9]). 
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To  characterize  the  external  or  active  patch  contributions,  it  is  typical  to  start  with  the 
assumption  that  the  strains  generated  by  a  patch  are  proportional  to  the  applied  voltage 
[3].  Since  differing  voltages  can  be  applied  to  the  outer  and  inner  patches  in  the  pair,  we 
will  differentiate  between  the  two  with  and  Vi2{t)  used  to  denote  the  voltages  to  the 

outer  and  inner  patches  in  the  pair,  respectively.  The  proportionality  constant  relating 
the  generated  strain  to  the  input  voltage  is  designated  by  d^i.  As  detailed  in  [3],  the  total 
external  moments  and  forces  generated  by  the  patches  are 


(M.U  = 


—  Epe  C?3i 


P^i  1—1/  h 


l^Xpe.{x,0) 


^  y.  as 

2  ^  3r)  V  2  ZR 


Vi 


i2 


(M.),,,  =  — 


<^310-2 


l^pe  ^hpQ 


XpeX^jmi-V,2] 


Epe  d^X 


1  ^pe  ^pe 


Xpei  (^5  ^)^pei 


^pe  H” 


^2 

2R 


^  Kl  +  (^hpe  — 


«2 

2R 


Vi 


42 


(2.4) 


(Ne) 


-E, 


pe 


pei 


V, 


^SlXpe;  (^)  [ftl  Vi2 


pe 


where  a2  =  {h/2  +  hpef  -  h^f4  and  Ug  =  {h/2  +  hp^f  -  h^/S.  When  substituted  into  (2.1), 
the  expressions  (2.4)  provide  the  input  from  the  patches  when  voltages  are  applied. 

Finally,  the  fixed-edge  boundary  conditions 


u 


V  =  w 


^  =  0  ,  x  =  0,£  (2.5) 

are  used  to  model  the  end  behavior  of  the  shell.  These  boundary  conditions  are  appropriate 
for  experimental  setups  in  which  heavy  endcaps  prevent  edge  movement.  Note  that  alternative 
boundary  conditions  such  as  simply  supported  or  “almost  fixed”  (see  [4])  can  be  employed  if 
edge  movement  is  suspected. 


2.2  Weak  Form  of  Modeling  Equations 

The  strong  form  (2.1)  of  the  modeling  equations  requires  first  and  second  derivatives  of  the 
moment  and  force  resultants.  As  noted  in  (2.2)  and  (2.4),  both  the  internal  and  external 
moment  and  force  resultants  are  discontinuous  due  to  the  piezoceramic  patches.  Hence  formal 
analysis  and  approximation  using  the  strong  form  of  the  modeling  equations  lead  to  difficulties 
due  to  differentiation  of  Dirac  distributions. 

To  alleviate  these  difficulties,  it  is  advantageous  to  consider  a  weak  form  of  the  modeling 
equations  which  can  be  derived  from  Hamilton’s  principle  (energy  considerations).  While 
equivalent  to  the  strong  form  under  suitable  smoothness  assumptions,  the  weak  form  provides 
a  more  natural  setting  for  analysis  and  approximation. 

The  state  variables  for  the  problem  in  second-order  form  are  taken  to  be  j/  =  (u,u,u))  in 
the  state  space  H  =  T^(ro)  x  T^(ro)  x  T^(ro).  For  the  fixed-edge  boundary  conditions  (2.5), 
the  space  of  test  functions  is  taken  to  be  F  =  Hq{To)  x  i/o(ro)  x  ifo(ro)  where 

^^o(ro)  =  {r?  e  if^ro)  h(O,0)  =  =  0} 

^^0  (ro)  =  {ve  H\To)  I  v{0, 0)  =  7^,(0,  e)  =  6)  =  0)  =  0}  . 
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For  $  =  (u,v,w)  and  $  =  (?/i,  772,  ^ys))  the  H  and  V  inner  products  are  taken  to  be 

/  phuTiid'y+  /  phvrj2d'y+  /  phw%d'y 
JVq  •'To  ‘'To 

and 


i=l  ^  ‘^pe 


/  \  .  dr]i 

{e^  +  ’^pe£e)-^  de 
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/  1  /1  \  ^^2 


"  2Epehpe 


Xpci 


''pe 


,  ,^772  .  I  .  dr]2 

{ee  + 


^7 
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r  f_^ 
Jr.  1 1  -  i' 


I  ,  , _  .  ^^7/3 

-{eg  +  uea;)p3  + 


/i^  ,  ,  ^d'^Vs  ^  d'^m 


12i?2 


dxdO 
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3(1  -  7/2J 


Xpei  (^7  ^) 


^^pe  _  \:rr—  ^  1  .. 

"1“  ^pe^xjV^  "I"  ^pe^6) 


■  j^2^i^e  +  ^pei^x) 


■  d-y 


where  ex^eg,ex9,  Hx,  Kg,T  are  defined  in  (2.3)  and  ^7  =  RdOdx.  The  dependence  of  the  inner 
product  on  the  Young’s  moduli  is  explicitly  included  in  the  definition  to  provide  a  notation 
for  defining  analogous  damping  expressions  later  in  this  work. 

The  weak  form  of  (2.1),  as  derived  in  [5]  from  energy  principles,  is  given  by 

L  {“£•'  *  "-5  ♦  “••w  -  -  «b'<.u£  j*  ■  • 


L 


Rph^V2  +  -  Rgn,  ^  )d-,  =  0 


de 


dx 


i=l 


de 


f  f ,d'^w_  „  dw_  d‘^r)3  1 


d^m  1  „  a^.,3 


de^ 


dxde 


—  RqnTI^  + 


j=l  L 


R{Mx) 


,  1 


p^i 


dx^  R 


+  -B(Me) 


pet 


d^Vs 

de^ 


>dj  =  0 


(2.6) 


6 


for  all  ^  =  (?7i,»?2,»73)  G  V.  A  comparison  between  (2.6)  and  (2.1)  illustrates  that  in  the 
weak  form,  derivatives  are  transferred  from  the  discontinuous  resultants  onto  suitably  smooth 
test  functions.  This  alleviates  the  difficulties  associated  with  the  discontinuities  and  reduces 
smoothness  requirements  on  approximate  solutions. 

2.3  Abstract  Formulation 

To  define  appropriate  sesquilinear  forms,  we  group  stiffness  components  separately  from  damp¬ 
ing  components.  To  this  end,  we  define  <7*-  :  V"  x  1/  — >  C  ,  i  =  1,2  by 

(T2($,’I')  =  ^(cD,cnpJ$,^)  +  /  nwfj^d') .  ^  ^ 

^  «/ro 

Note  that  ((c^i,  ^)y  differs  from  {{E,  Epe)^,'9)y  only  in  that  Young’s  moduli  are 

replaced  by  Kelvin- Voigt  damping  coefficients.  It  can  be  directly  verified  that  the  stiffness 


form  (Ti  satisfies 

(HI) 

|<Ti($,ff')|  <  Ci|#|y|^|y  ,  for  some  Cl  G  IR 

(Bounded) 

(H2) 

Reai(#,$)  >  C2|$|y  ,  for  some  C2  >  0 

(y-Elliptic) 

(H3) 

ai($,^)  =  C7i(^,$) 

(Symmetric) 

for  all  €  V. 

Moreover,  the  damping  term  a2  satisfies 

(H4) 

|cr2($,  ff')!  <  Csl^lyl^'jy  ,  for  some  C3  G  IR 

(Bounded) 

(H5) 

Recr2($,$)  >  C4|$|y  ,  for  some  C4  >  0 

(y-Elliptic) . 

Remark  1.  The  symmetry  of  ai  is  dependent  upon  the  choice  of  shell  model  and  ultimately 
reflects  the  Maxwell-Betti  reciprocity  theorem.  While  the  Donnell-Mushtari  model  yields  a 
symmetric  sesquilinear  form  (Ti,  other  models  such  as  the  Timoshenko  shell  model  will  not 
yield  a  symmetric  form. 

To  represent  control  contributions,  let  U  =  ]R^  denote  the  Hilbert  space  of  control  inputs 
and  define  B  ^  C{U,  V*)  by 

for  ^  e  y.  Here  (•,  ■)y,  y  denotes  the  usual  duality  product.  Finally,  with  the  definition 
g  =  {Ij ph)[qx,q6.,qn],  we  can  write  the  weak  form  (2.6)  in  the  abstract  variational  form 

^)  +  ^)  =  {Bu{t)  -|-  g{t),  •  (2.8) 
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To  pose  the  problem  in  a  first-order  form  amenable  for  control  applications,  we  define  the 
product  spaces  'H  =  V  x  H  and  V  =  V  xV  with  the  norms 

+  \^2\h 

1(01,  <^2) Iv  =  I0ily  +  \4>2\v  • 


The  state  is  taken  to  be  z{t)  =  {y{t),y{t))  e  Finally,  the  product  space  forcing  terms  are 
formulated  as  ,  , 


9{i)  = 


0 


Bu{t)  — 


0 

Bu{t) 


The  weak  form  (2.8)  can  then  be  rewritten  as 

{i((),A)y.  y  +<7(2((),A)  =  {Btl(l)  +  j((),A)y.  y  foP  A  €  V 
2(0)  =  20  =  (vo,Vi) 


(2.9) 


(2.10) 


where  cr ;  V  x  V  — >•  €  is  given  by 


(t(0,0)  =  -  (02,0l)v  +  '^l(01’02)  +  0-2(02,02) 


for  0  =  (0i,02),0  =  (01,02)  ^  V.  As  proven  in  [5,  page  109],  <t  is  V  continuous  and  for 
A  >  0,  (t(-,-)  +  V-elliptic.  From  the  continuity  of  cr,  it  follows  that  one  can  define 

an  operator  A  G  £(V,  V*)  by  <j(T,A)  =  (.4T,  A)v»,v  • 

To  obtain  a  strong  form  of  the  first-order  system  which  is  appropriate  for  control  purposes, 
consider  the  system  operator 


dom^  =  {(01,  02)  ^  ^102  ^  F,  Ai0i  -f  ^202  ^  ff} 

-Ai  -A2 


(2.11) 


with  Ai,A2  €  C{V^V*)  defined  by 

(Aj0i,  02)y»^y  —  <^i(01,02)  ,  i  =  1,2  . 

It  should  be  notated  that  A  is  the  negative  of  the  restriction  to  dom.4  of  ,4  €  £(V,V*)  so 
that  (7(T,A)  =  (-.4T,A)„  for  T  G  dom^A  G  V. 

A  strong  form  of  the  abstract  system  model  is  then  given  by 

z{t)  =  Az{t)  -k  Bu{t)  +  g{t)  in  V*  =  y  X  y*  ^ 

•s(O)  =  zo . 

The  rigorous  equivalence  of  solutions  is  established  through  the  following  theorems. 


Theorem  1.  Under  Hypotheses  (H1)-(H5)  on  <7i  and  <72,  A  generates  an  analytic  semigroup 
T{t)  on  V,'H  and  V*.  In  terms  of  this  semigroup,  the  representation 

z{t)^T{t)zo+  [  T{t  -  s)[Bu{s)  +  g{s)]ds  (2.13) 

Jo 

defines  a  mild  solution  to  (2.12)  for  zq  G  V*  and  Bu  +  g  e  i/^((0,  T);  V*).  Furthermore,  this 
semigroup  is  (uniformly)  exponentially  stable  on  V,7i  and  V*. 
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Theorem  2.  Let  Zgg  denote  the  semigroup  solution  to  (2.12)  given  by  (2.13)  and  let  v^ar  denote 
the  weak  solution  to  (2.8).  Under  hypotheses  (H1)-(H5),  it  follows  that  Zsg{zo,  T)  =  Zyar{zoi  T) 
for  zo  eH  and  !F  =  Bu  +  g  £  L^((0,  T);  V*). 

Following  the  convention  of  [14],  we  will  use  the  same  notation  for  the  semigroups  defined 
onV,'H  and  V*  since  each  semigroup  is  an  extension  or  re^riction  of  the  others.  Note  that 
dom.4  defined  in  (2.11)  is  actually  dom^^.^,  the  domain  of  ^  as  a  generator  of  T{t)  in  As 
detailed  in  Lemma  3.6.1  and  Theorem  3.6.1  of  [14]  (see  also  Section  IV.6  of  [12]  and  Chapter  2, 
Theorem  5.2  of  [11]),  the  property  that  A  generates  an  analytic  semigroup  on  V,'H  and  V* 
results  from  the  continuity  and  V-ellipticity  of  a.  The  exponential  stability  of  T{t)  on  V. 
for  second-order  systems  with  strong  damping  is  demonstrated  in  [1]  while  the  exponential 
stability  of  T{t)  on  V  and  V*  in  this  setting  is  proven  in  Lemma  3.3  of  [2].  Finally,  Theorem  2 
is  a  reformulation  of  Theorem  4.14  of  [5]  and  details  can  be  found  therein. 

3  LQR  Control  Problem 

In  the  last  section,  the  PDF  system  modeling  the  dynamics  of  the  thin  shell  with  surface- 
mounted  piezoceramic  actuators  was  written  in  the  abstract  first-order  form 

z{t)  =  Az{t)  4-  Bu(i)  -|-  g(t) 

z(0)  =  zo 

in  V*.  In  this  section,  LQR  control  results  for  both  the  original  infinite  dimensional  problem 
and  approximating  finite  dimensional  problems  will  be  discussed.  Two  cases  will  be  consid¬ 
ered,  namely  when  g  =  0  and  g  is  periodic  in  time.  In  both  cases,  it  is  assumed  that  state 
observations  in  an  observation  space  Y  have  the  form 

Zob{t)  =  Cz{t)  (3.1) 

where  C  €  C{'H,  T)  is  bounded.  The  assumption  that  C  is  bounded  is  made  to  simplify  the 
exposition  and  the  reader  is  referred  to  [2]  for  arguments  pertaining  to  the  case  of  unbounded 
observation  operators. 


3.1  No  Exogenous  Input 

For  the  case  in  which  ^  =  0,  the  infinite  horizon  problem  concerns  the  determination  of  a 
control  u  which  minimizes  the  quadratic  cost  functional 

J{u,zo)  =  j^  {\Cz{t)\^  +  \TZ^/\{t)\'lj'jdt  (3.2) 

subject  to 

z{t)  =  Az{t)  Y  Bu{t) 

^(0)  =  Zo . 

The  positive,  self-adjoint  operator  TZ  is  used  to  weight  various  components  of  the  control. 
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As  detailed  in  [2,  5],  if  {A,B)  is  stabilizable  and  {A,C)  is  detectable,  then  the  optimal 
control  minimizing  (3.2)  is  given  by 


u{t)  = 

where  If  solves  the  algebraic  Riccati  equation 

{A*Il  +  IIA-  +  C*C)z  =  0  for  all  z  G  V 

and  z{t)  =  S{t)zo.  Here  S{t)  is  the  closed  loop  semigroup  generated  by  ,4  —  B7i~^B*U. 

For  implementation  purposes,  it  is  necessary  to  define  an  approximate  system  and  controls, 
and  determine  convergence  criteria  for  these  approximate  controls  when  fed  back  into  the 
infinite  dimensional  system.  The  approximations  are  considered  in  a  Galerkin  framework 
with  trajectories  evolving  in  the  finite  dimensional  subspaces  C  V  C  'H.  It  is  assumed 
that  the  approximation  method  satisfies  the  standard  convergence  conditions 

(HIN)  For  any  2  G  V,  there  exists  a  sequence  G  such  that  \z  —  z^\v  — >  0  as  A  — >  00. 

The  finite  dimensional  operators  and  approximating  system  are  then  determined  as  follows. 
The  operator  A^  :  which  approximates  A  is  defined  by  restricting  cr  to  x 

this  yields 

{-A^T,  a)^  =  <t(T,  A)  for  all  T,  A  G  .  (3.3) 

For  each  A,  the  Cq  semigroup  on  which  is  generated  by  ,4^  is  denoted  by  it).  The 
control  operator  is  approximated  by  G  £(17,  V^)  given  by 

(B^u,a)  =(u,H*A)^  foralluGl7,AG  (3.4) 

\  f  'hi 

while  denotes  the  restriction  of  the  observation  operator  C  to  V^.  Finally,  we  let 
denote  the  usual  orthogonal  projection  of  Ti  onto  which  by  definition  satisfies 

(i)  P^T  G  for  T  G  -H 

(ii)  (P^T  -  T,  A)^  =  0  for  all  A  G  . 

This  projection  can  be  extended  to  P^  G  £(V*,  V'^)  by  replacing  the  'H-inner  product  (T,  A)^ 
by  the  duality  product  (T,  A)y,  y  and  considering  T  G  V*. 

The  approximate  problem  corresponding  to  (2.10)  with  g  =  0  can  then  be  formulated  as 

7  a)^  +  <7(z~(0,A)  =  {B''u(<),a)^  tor  all  A  €  V" 

^JV(o)  ^  _ 

This  has  the  solution 

z^{t)  =  r^{t)P’^zo  +  /'V^(f  -  s)P^B^u{s)ds  . 

Jo 
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The  following  theorems  taken  from  [2,  5]  can  be  used  to  establish  the  convergence  of  the 
approximate  gains  to  their  infinite  dimensional  counterparts  for  certain  classes  of  shell  models 
(see  specifically  Theorem  7.10  and  Lemma  7.13  of  [5]). 

Theorem  3.  Assume  that  the  injection  i  :  V  H  is  compact.  Moreover,  suppose  that  the 
damping  sesquilinear  form  can  be  decomposed  as  <T2  =  Sai  +  a2,  for  some  ^  >  0,  where  the 
continuous  sesquilinear  form  da  satisfies  for  some  A  G  IR 

Re  (T2(^,  0)  ^  2 1*^1^  4>  G  V  . 

Finally,  suppose  that  the  operator  A2,  where  A2  G  £(V,  F*)  is  defined  by  ^  = 

(72(0,  ?7),  is  compact  on  V. 

If  for  some  a;  G  IR  and  M  >  1,  T{t)  satisfies 

\T{R)\c{h)  <  t  >  0  , 

then  for  any  £  >  0  there  exists  an  integer  such  that  for  N  >  N^, 

|r^(t)P^U(«)  <  ,  t>0 

for  some  constant  M  >  0  independent  of  N. 


Theorem  4.  Assume  that  the  injection  i  :  V  H  is  compact.  Let  the  sesquilinear  form  cr 
associated  with  the  first-order  system  (2.10)  be  continuous  and  V-elliptic.  Assume  that  the 
operators  A^B^C  of  (2.11),  (2.9),  (3.1),  respectively,  satisfy:  {A,B)  is  stabilizable  and  (A.,C) 
is  detectable  where  B  G  C{U,V*)  is  unbounded  and  C  G  €{71,  Y)  is  bounded.  Consider  an 
approximation  method  which  satisfies  (HIN).  Finally,  suppose  that  for  fixed  No  and  N  >  No, 
the  pair  [A^,B^')  is  uniformly  stabilizable  and  {A^,C^)  is  uniformly  detectable. 

Then  for  N  sufficiently  large,  there  exists  a  unique  nonnegative  self-adjoint  solution  11^  G 
jC{V*,  V)  to  the  approximate  algebraic  Riccati  equation 

-f  u^A^  -  =  o 

in  V^.  There  also  exist  constants  M3  >  1  and  (^3  >  0  independent  of  N  such  that  S^it)  = 


{%,  <  ,  f>0. 


Moreover,  the  convergence  of  the  Riccati  and  control  operators 

-4-  in  V  for  every  2^  G  V* 


pN 


CCKfJ) 


as  N  -¥  00,  is  obtained. 
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Example  1.  We  consider  in  this  example  a  shell  with  constant  parameters  p,  E,  u,  c©.  Such  a 
case  would  arise  if  modeling  a  homogeneous  shell  or  a  shell  in  which  the  variance  of  material 
properties  across  regions  with  actuators  is  negligible.  The  sesquilinear  forms  for  this  model 
are  specified  in  (2.7).  Due  to  the  constant  coefficients,  <T2  can  be  written  as  <T2  =  Sai  +  a-2 
where  S  =  ^  and  (J2 ($,'!')  =  p  /pg  wrjzd'^.  It  follows  immediately  that 

Red2{(t>,4>)  =  1^  ^ 

for  all  (j)  e  V.  The  boundedness  of  the  operator  A2  generated  by  (T2  follows  directly  from 
the  boundedness  of  0-2.  Furthermore,  it  is  noted  that  e  C{V*,  V)  can  be  written  as  an 
operator  on  V  -4  V  by  A^^  =  Ai^i*i  where  the  injections  ij  V  H,i*  :  H  M-  V*  are 
compact.  Thus  A^^  is  compact  on  V  which  implies  that  A7M2  is  compact  on  V  since  it  is 
formed  from  the  product  of  compact  and  bounded  linear  operators.  Finally,  the  exponential 
stability  of  T{t),  the  stabilizability  of  {A,  B)  and  the  detectability  of  {A,  B)  are  guaranteed  by 
Theorem  1.  The  hypotheses  of  Theorem  3  are  then  satisfied  for  this  system  and  one  obtains 
uniform  bounds  on  the  approximating  semigroups.  The  convergence  of  the  Riccati  and  control 
operators  is  then  obtained  from  Theorem  4. 

3.2  Periodic  Exogenous  Input 

A  reasonable  assumption  in  many  mechanical  systems  is  that  g  is  periodic  in  time  with  period 
r.  The  system  to  be  controlled  in  this  case  is 

z{t)  =  Az{t)  A  Bu{t)  +  g{t) 

(3.5) 

^(0)  =  z{t) 

and  an  appropriate  quadratic  functional  to  be  minimized  is 

L  JO 

Note  that  the  periodic  exogenous  term  g  can  be  used  to  model  inputs  such  as  noise  gener¬ 
ated  by  rotating  engine  components  (e.g.,  propellers  or  turbines)  or  periodic  electromagnetic 
disturbances. 

To  guarantee  the  existence  of  a  unique  Riccati  solution  and  control  for  the  system  (3.5), 
it  is  assumed  that  {A,B)  is  stabilizable  and  {A,C)  is  detectable.  Furthermore,  it  is  assumed 
that  g  e  L^(0,r;'H)  and  that  B  is  bounded.  Under  these  conditions,  it  is  verified  in  [6]  that 
the  Riccati  equation 

A*U  +  UA  +  UBn-^B*Il  +  CX  =  0 

has  a  unique  solution.  Furthermore,  if  r  denotes  the  r-periodic  solution  of  the  adjoint  or 
tracking  equation 

r{t)  =  -[A-Bn-X*TlYr{t)  +  Ilg{t) 
r(0)  =  r(T) 
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and  2:  is  the  closed  loop  solution  of 


i(t)  =  IA-  eK-‘e*n]s(()  -  sn-'B'r(t)  +  g(t) 

z(0)  =  f(T), 

then  the  optimal  control  is  given  by 

u{t)  =  -  r{t)] .  (3.6) 

The  LQR  theory  for  this  case  is  less  complete  than  that  for  systems  with  no  exogenous 
input  and  is  currently  limited  to  bounded  control  inputs  B.  The  synthesis  of  the  theory  for 
unbounded  input  operators  and  periodic  exogenous  forces  is  currently  under  investigation. 
The  effectiveness  of  the  method  is  illustrated  in  the  final  example  of  this  work. 


4  Approximation  Method 

A  Galerkin  method  was  used  to  approximate  the  solutions  u,v,w  to  the  system  (2.6),  or 
equivalently,  (2.10).  The  approximating  subspaces  were  taken  of  the  form  =  span{H„^}x 
span{0„^}  X  span{Bu,^}  where  Bu^,Bv^,Bwi,  denote  bases  for  the  u,v  and  w  displacements, 
respectively.  To  exploit  the  tensor  nature  of  the  shell  domain  To  and  periodicity  in  6,  the  bases 
were  constructed  with  Fourier  components  in  9  and  cubic  splines  in  x  (see  [7]  for  details).  The 
approximate  displacements  were  then  given  by  the  expansions 

u^{t,  6,x)  =  Y^  Uk{t)Buk{9,  x) 

A;=l 

v^{t,e,x)  =  ^Vk{t)B,,{9,x) 

k=i 

w^{t,9,x)  =  J2'^k{t)Bn,^{6,x) . 

k=l 

To  obtain  a  finite  dimensional  system  with  matrices  corresponding  to  the  finite  dimensional 
operators  in  (3.3)  and  (3.4),  the  sesquilinear  forms  ai  and  (72  were  restricted  to  V-^.  This  yields 
the  matrix  system 


0  ■ 

’  9^  (t)  ' 

0 

1 

■  9^  it)  ' 

0 

[^t)]  + 

■  0  ■ 

0 

. 

_  {t) 

1 

1 

0^  1 

9^{t) 

+ 

.  9^(()  . 

0  ■ 

'  i9^(0)  ■ 

■  yf  ' 

0 

??^(0) 

.  y2  . 

where  ??^(f)  _  [^^(t),  •  • . ,  Uj^^,vi(t),  •  •  ■ ,  contains  the  M  =  Mu+Mv+Myj 

generalized  Fourier  coefficients.  The  s  patch  inputs  are  contained  in  u{t)  —  [ui(f),  •  •  • ,  Us(f)]^. 
The  reader  is  referred  to  [7]  for  details  concerning  the  construction  of  the  mass,  stiffness  and 
damping  matrices  M'^,  ,  Kf^,  the  inputs  and  the  initial  conditions  yf  . 
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(4.1) 


Multiplication  by  the  inverted  mass  matrix  yields  the  Cauchy  equation 

z^{t)  =  A^z^{t)  +  B^u{t)  +  g’^it) 

AO)  = 

where  z^  €  This  system  forms  the  constraint  equations  used  in  the 

finite  dimensional  LQR  theory  discussed  in  Section  3. 

5  Numerical  Example 

We  consider  here  an  exogenous  force  g  which  is  periodic  in  time  with  period  r  =  IOOOtt 
(500  Hz).  The  distribution  of  the  force  was  taken  to  be  binormal  in  the  transverse  and 
longitudinal  directions  and  centered  at  {x,0)  =  (^/2,0)  and  {x,9)  =  (^/2,  tt)  as  depicted  in 
Figure  2.  The  magnitude  of  the  transverse  component  qn  was  one  hundred  times  that  of  the 
longitudinal  component  qx  so  as  to  model  an  input  consisting  primarily  of  acoustic  sources 
located  adjacent  to  {i/2,0)  and  {i/2,7r). 

Six  pairs  of  piezoceramic  patches  of  length  1  cm  and  radial  measure  7r/3  were  employed 
as  actuators.  The  locations  and  material  properties  of  the  patches  along  with  the  dimensions 
and  physical  parameters  for  the  shell  are  summarized  in  Table  1. 

To  accommodate  the  periodic  exogenous  force  g,  control  inputs  to  the  twelve  patches  were 
computed  using  the  feedback  law  (3.6).  Note  that  in  this  formulation,  independent  voltages 
are  determined  for  the  individual  patches.  This  provides  the  capability  of  generating  both 
inplane  forces  and  bending  moments  in  the  regions  covered  by  the  patches  so  that  longitudinal, 
circumferential  and  transverse  vibrations  can  be  controlled. 

Time  histories  of  the  uncontrolled  and  controlled  shell  displacements  at  the  point  pi  = 
(3^/4,;r/32),  depicted  in  Figure  2,  are  plotted  in  Figure  3.  The  open  loop  trajectories  exhibit 
both  a  transient  response  settling  into  steady  state  and  a  beat  phenomenon  due  to  the  close 
proximity  of  the  driving  frequency  and  natural  frequencies  for  the  shell.  At  this  observation 


Figure  2.  Distribution  of  normal  forcing  function  at  0  =  0  and  $  =  w.  Observation  lines 
Li  =  {(2;, 0)10  <x<i,6^  7r/6},  L2  =  {(a;,0)|.T  =  3^/4, 0  <  0  <  27r}  and  observation  point 
Pi  =  (3^/4,7r/32). 
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point,  all  three  displacement  components  are  reduced  by  more  than  90%  when  controlling 
voltages  are  fed  back  to  the  patches. 

To  illustrate  the  spatial  attenuation  due  to  the  feedback  of  voltages  to  the  patches,  root 
mean  square  (rms)  plots  of  the  uncontrolled  and  controlled  trajectories  along  the  axial  line 
Li  and  circumferential  line  L2  (see  Figure  2)  are  plotted  in  Figure  4  and  5,  respectively. 
For  the  open  loop  case,  these  plots  illustrate  a  standing  wave  in  all  three  components  of  the 
displacement  (the  slight  asymmetry  in  the  axial  plots  is  due  to  the  longitudinal  input  q^). 
The  figures  also  demonstrate  significant  reductions  in  all  three  displacement  levels,  even  in 
regions  not  covered  by  patches.  This  further  illustrates  the  effectiveness  through  which  the 
model-based  control  law  can  be  used  to  attenuate  shell  vibrations. 


Time  (seconds) 

(c) 

Figure  3.  Uncontrolled  and  controlled  shell  displacements  at  the  point  pi  =  (3£/4, 7r/32); 

(a)  longitudinal  u,  (b)  circumferential  u,  (c)  transverse  w  displacements; - (uncontrolled), 

—  (controlled). 
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Figure  4.  Root  mean  square  (rms)  displacements  along  the  axial  line  Li] - (uncontrolled) 

—  (controlled). 


Figure  5.  Root  mean  square  (rms)  displacements  along  the  circumferential  line  L2] 
- (uncontrolled),  —  (controlled). 


Dimensions 

Parameters 

Shell 

h  =  . 00127  m 

R  =  Am 

i  =  Im 

p  =  2700  kg/m^ 
E^7.1x  10^°  iV/m2 

C£)  =  1.47  X  10®  Nms 
ly  =  .33 

p  =  58.97 Af5/m2 

Patches 

hpe  =  .0001778  m 

Centers  (xj):  (.25, 0),  (.5, 0),  (.75, 0) 
(.25,7r),(.5,7r),(.75,7r) 
Dimensions:  x  :  0.1  cm,  0  :  n/3 

Ppe  =  7000  kg /wA 

Epe  =  6.3  X  10^°  N/m^ 

—  1-7  X  10®  Nms 
t^pe  —  *31 

4i  =  190  X  lO-^^mlV 

Table  1.  Dimensions  and  physical  parameters  for  the  shell  and  patches. 


6  Conclusions 

A  model-based  LQR  method  for  controlling  shell  vibrations  has  been  presented  here.  While 
developed  in  the  context  of  a  modified  Donnell-Mushtari  cylindrical  shell  model,  the  method 
is  quite  general  and  can  be  directly  extended  to  other  models  and  geometries.  Under  the 
assumption  of  strong  or  Kelvin- Voigt  damping  (a  reasonable  and  typical  assumption  for  many 
shell  materials),  model  well-posedness  and  convergence  of  control  gains  is  obtained  using 
analytic  semigroup  theory  combined  with  LQR  results  for  unbounded  input  operators. 

The  Galerkin  method  used  to  approximate  the  system  dynamics  utilizes  bases  constructed 
from  tensored  Fourier  polynomials  and  modified  cubic  splines.  As  discussed  in  [7],  care  must 
be  taken  when  developing  methods  for  approximating  shell  dynamics  so  as  to  avoid  shear  or 
membrane  locking.  One  manifestation  of  locking  is  the  existence  of  model  dynamics  which 
are  incorrectly  approximated  by  the  numerical  method.  The  use  of  a  numerical  method  which 
exhibits  locking  can  lead  to  a  loss  of  control  authority  and  potential  controller  destabilization 
if  the  approximations  are  sufficiently  inaccurate.  Further  details  regarding  issues  concerning 
the  approximation  of  shell  dynamics  and  convergence  properties  of  the  numerical  method  can 
be  found  in  [7]. 

The  numerical  example  demonstrates  the  effectiveness  of  the  model-based  control  method 
for  attenuating  all  three  components  of  the  shell  displacement  in  the  presence  of  both  transient 
and  steady  state  dynamics.  Furthermore,  by  modeling  the  global  shell  dynamics  and  patch 
interactions  through  coupled  PDF  and  constructing  the  control  law  in  terms  of  these  PDF, 
significant  reductions  in  displacement  levels  throughout  the  shell  are  obtained,  even  in  regions 
devoid  of  patches.  Numerical  implementation  of  the  LQR  method  in  this  manner  provides  a 
first  step  toward  the  development  of  model-based  state  estimators  and  compensators  which 
can  be  experimentally  implemented  in  shell  applications. 
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